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ABSTRACT 

A theoretical analysis is presented to study the 
extinction characteristics of a diffusion flame near 
the leading edge of a thin fuel plate in slow, forced 
convective flows in a microgravity environment. The 
mathematical model includes two-dimensional Navier- 
Stokes momentum, energy and species equations with 
one-step overall chemical reaction using second-order 
finite rate Arrhenius kinetics. Radiant heat loss on 
the fuel plate is applied in the model as it is the 
dominant mechanism for flame extinguishment in the 
small convective flow regime. A parametric study 
based on the variation of convective flow velocity, 
which varies the Damkohler number (Da), and the sur- 
face radiant heat loss parameter (S) simultaneously, 
is given. An extinction limit is found in the regime 
of slow convective flow when the rate of radiant heat 
loss from fuel surface outweighs the rate of heat 
generation due to combustion. The transition from 
existent envelope flame to extinguishment consists of 
gradual flame contraction in the opposed flow direc- 
tion together with flame temperature reduction as the 
convective flow velocity decreases continuously until 
the extinction limit is reached. A case of flame 
structure subjected to surface radiant heat loss is 
also presented and discussed. 

NOMENCLATURE 


E activation energy for gas phase reaction, E/RT„ 
Es activation energy for fuel pyrolysis, Es/RT„ 
f stoichiometric oxidizer/fuel mass ratio 

k* thermal conductivity at T* [W • M^k -1 ] 

L latent heat, L/cpT„ 

Le Lewis number, a/D 

m blowing rate on the fuel plate, m/p*u„ 

Pr Prandtl number. Via 

p pressure, (p - p rc )/p*u£ 

p rc absolute pressure at arbitrary reference 

point [Pa] 

Q heat of reaction per unit mass of fuel, Q/CpT„ 

Q rate of heat generation per unit volume, Quf 

R universal gas constant [kg • m 2 • S -2 • kg - 
mole - '*- • k -1 ] 

Re Reynolds number, p*u„5/u* 


As pre-exponential factor for fuel pyrolysis. As/p* 

B dimensional frequency factor for gas phase 
reaction, [m 2 • kg - *- • S - *-] 

cp average specific heat [kj • kg -2 • k -2 ] 

D dimensional species diffusivity [m 2 • S -2 ] 

Da Demkohler number, (k*/'cpu 2 )B 
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r ratio of reference temperature to ambient 

temperature, T*/T„ 

S surface radiant heat loss parameter, coT 2 / e *CpU 

T* reference temperature [K] 

T temperature, T/T„ 

T s pyrolysis temperature, Ts/T„ 

T w fuel surface temperature, T w /T m 

T„ ambient temperature [k] 



u velocity parallel to fuel surface, u/u„ 

u„ free stream velocity, [m • S - *] 

v velocity normal to fuel surface, 'v/u„ 

x distance along fuel surface, "x/a 

x max location of upstream boundary 

Xffljp location of downstream boundary 

Yp mass fraction of fuel, pp/'p 

Yp„ ambient fuel mass fraction 

Yg mass fraction of oxidizer, og /'p 

Yg„ ambient oxidizer mass fraction 

y distance normal to fuel surface, y Is 

Ymax location of free stream boundary 

a* thermal diffusivity at T*, [m^ • S - *] 

6 thermal distance, a*/u„, [m] 

e total hemispherical emissivity of the fuel 
surface 

p* dynamic viscosity at T*, [kg * nr 1 • S -1 ] 
u dynamic viscosity, H/v* 

p* density at T*, [kg • m~3] 

p density, p/p* 

o Stef an-Bol tzman constant [W • • K -/ *] 

up nondimensional fuel reactivity, -Dap 2 YpYgEXP(-E/T) 

Overhead 

dimensional quantities 
Superscript 
* reference state 

Subscript 
F fuel 

0 oxidizer 

w surface of fuel plate 

ambient 
max maximum 

min minimum 

INTRODUCTION 

The study is concerned with the extinction mech- 
anism of diffusion flame adjacent to a thin fuel plate 
in slow convective flow as shown in Fig. 1. The major 
application for this class of problem is fire safety 
in manned spacecraft and space station where the 


gravity effect is negligible. Under these circum- 
stances, natural convection is neglected. But a 
forced flow can exist due to the air ventilation in 
the living quarters of the astronauts. This convec- 
tive velocity can be small, less than the natural- 
convective velocity around the flame in earth 
environment which is of the order of 15 cm/s. With 
the increasing concern of fire safety in spacecraft 
and accelerated experimental efforts to test the 
flammability characteristics of solid materials in 
microgravity at NASA Lewis Research Center, this 
effort describes the theoretical investigation coor- 
dinated with experimental task. 

The present analysis is an extension of earlier 
work U ) . In that work, the blowoff phenomena for a 
two-dimensional diffusion flame around the leading 
edge of a thin fuel plate were investigated. The 
primary parameter varied in the computation is the 
Damkohler number. Da, defind as (p*a*/u£)B where 
p* and a* are density and thermal diffusivity of 
the gas in reference state respectively, u„ is the 
convective flow velocity and B is the frequency 
factor of chemical reaction. Lowering Da corre- 
sponds to increasing u„ for a given fuel. The blow- 
off phenomena predicted by the previous theoretical 
model (1_) is what we observe in ordinary experience. 
However, the velocity range of interest in micro- 
ravity application is when u„ becomes small. In 
ef. 1, only one nondimensional parameter contains 
u„, i.e., the Damkohler number. When u„ * 0, Da ♦ «. 
which according to (l_).is moving away from blowoff or 
extinction. At the same time, it is also found that 
dimensionally the fuel burning rate is very small and 
the flame stand-off distance is large. In other 
words, in the limit of small convective flow velocity, 
a kinetically-strong (large Da), large size but low 
power (small burning rate) flame exists. This con- 
clusion is a contradiction to the experimental obser- 
vations (2,3) which show lower flame temperature and 
flame spreading rate in zero-g condition. It implies 
that certain key elements may be missing in the model 
for small convective flow regime. The missing link 
appears to be the radiant heat loss which becomes 
important for a flame with small heat release rate 
due to reduced oxidizer supply rate. In this case, 
heat loss owing to radiation can be the dominant 
mechanism for flame extinction. 

Bonne (4), to the author's best knowledge, is 
the first one to investigate the influence of radia- 
tive heat transfer on diffusion .flame under zero-g 
environment. From both experimental work and theoret- 
ical analysis by using a flat diffusion flame, he 
concludes that the radiant heat loss results in the 
reduction of flame temperature which leads to extin- 
guishment. Recently, T'ien (5) imposed a radiative 
heat loss term into the energy balance on the con- 
densed fuel surface for a stagnation-point diffusion 
flame. Without the effect of gravity, he identifies 
a new extinction branch at the regime of sufficiently 
small flame stretch rate and also finds the radiative 
contribution to extinction is through the reduction 
of flame temperature. 

In the other studies the primary purposes may not 
be microgravity combustion research, but they can give 
the important indication of the role of radiation in 
contributing to flame extinction. Sibulkin, Kulkarni 
and Annamalai (6) analyze the radiative effects on the 
laminar, natural-convective burning of a thermally 
thick, vertical fuel surface by a use of local simi- 
larity approximation. The results show that gas phase 
radiation has a negligible effect on burning rate, 
however, heat losses due to surface radiation are a 
sufficient condition for extinction even in the limit 
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of infinite-fast chemical reaction. In Ref. 7, 

Sohrab, Linan and Williams applied asymptotic tech- 
nique to investigate the diffusion flame with gaseous 
radiant heat loss in a stagnation-point boundary layer 
geometry. Another work by Sohrab and Williams (§_) 
studied the diffusion flames adjacent to flat, con- 
densed fuels with surface radiation heat loss. In 
both papers, the authors find that the diffusion 
flames in the stagnation-point flow can be extin- 
guished by providing additional radiant losses from 
the gas phase or on the fuel surface which will 
decrease the available enthalpy of combustion of the 
fuel, thereby lowering the flame temperature. 

The specific problem to be treated in this study 
is to apply the radiative heat loss on the fuel sur- 
face to the combustion model of Ref. 1. With the 
additional heat loss term in the energy balance, the 
burning rate of the fuel is expressed by an Arrhenius- 
type pyrolysis law. Also it is assumed that the flow 
is steady, laminar and two-dimensional, and no burnout 
of fuel sample is included. Radiative heat transfer 
in the gas phase is not considered here. 


Oxidizer species equation: 



( 6 ) 


where 

= -Dap 2 Y p Y 0 EXP (-E/T) (7) 

is the nondimensional fuel reaction rate for which the 
definition of Da can be found later in Eq. (16). 

Equation state: 

pT = y- = r (8) 

where r is the ratio of reference temperature to the 
ambient one (see Ref. 1). 


MATHEMATICAL MODEL 

Except for the modified boundary conditions, such 
as the radiant heat loss in the energy balance and the 
pyrolysis expression applied on the fuel surface, the 
present combustion model is the same as that in Ref. 1, 
The governing system is elliptic in nature which 
includes conservation equations for mass, momentum, 
energy and species in differential form. An equation 
of state and the variation of viscosity with tempera- 
ture are used as the linkages between properties. An 
assumed chemical reaction rate is also included in 
this system. Finally, a set of known boundary condi- 
tions is specified to complete the mathematical 
description of the problem. 

The nondimensional governing equations and bound- 
ary conditions, with the coordination defined' in 
Fig. 1, are as follows. 


The viscosity variation with temperature is 
taken to be 
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Boundary conditions: 
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Momentum equations: 
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Energy equation: 

ou 17 * c,v fy = HePF 


[h (■ S) * wj 


+ Qwp 


Fuel species equation: 
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KePr'fy|w = + S(Ts 4 - 1) 

where As and Es are a pre-exponential factor and 
activation energy for surface pyrolysis respectively. 
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S is the surface radiant heat loss parameter; 
defined as 

T 3 
eo I 

S = — — — - (15) 

p*C u 

P eo 

where e is the total hemispherical emissivity of the 
fuel surface and o is the Stefan-Boltzman constant. 

From above definition, S is a function of u„ 
if c /= 0. Another parameter in this study which also 
contains u„ is the Oamkohler number, shown in 
Eq. (7), expressed as 



For a given fuel under specified ambient condi- 
tion, changing the convective flow velocity (u„) will 
change both Da and S (see Eqs. (15) and (16)). 
Eventually, large Da implies a kinetical ly-strong 
flame in the gas phase. On the other hand, large S 
implies large heat loss to the fuel plate. So the 
resultant flame structure apparently depends on the 
competition between Da and S. 

One important feature of the boundary conditions 
at the fuel surface needs to be mentioned. Because of 
the use of Arrhenius pyrolysis law for the fuel, the 
temperature distribution along the wall is no longer ■ 
prescribed as in Ref. 1, but provided as a result of 
the calculation. 

The normalization procedure for these equations 
is described in Ref. 1. The problem is solved numer- 
ically. The computer code using the SIMPLE procedure 
(9) is developed and modified by the author. The 
algorithm is vectorized in order to utilize the per- 
formance of the supercomputer CRAY-l/S at Lewis. 

RESULTS and DISCUSSIONS 

There are 13 parameters in this study which are 
derived from the normalization procedure through the 
governing system. Their composite groups of physical 
properties and quantities are shown in Table I. The 
values of these properties are the same as that used 
in Ref. 1. As to the additional properties, such as 
As and Es, their values are adopted from Ref. 5. 

Only the effect of a variable convective flow velocity 
(tu), which changes Da and S simultaneously, is 
investigated with and without radiant heat loss ( e = 1 
and e = 0). According to the variation of the flow 
velocity in the regime of interest (<20 cm/s), the 
computation is categorized into 8 cases which are 
1 isted in Table II. 

For c = 0 which implies no surface radiant heat 
loss is involved, the results confirm that the simi- 
larity relationship still exists at large Da as 
found in Ref. 1, even though the surface boundary 
conditions for both studies are different. 

Figure 2(a) shows the temperature contour distribu- 
tions for Da ranging from 2.81x10® to 7.20x10®. 

For this range of Da envelope flames exist in the 
gas phase ahead of the leading edge (x = y = 0) of 
the fuel plate. Except for a very narrow region of 
high temperature, for example T > 6.9 (which are not 
shown in this figure), the temperature distributions 
are found to be nearly invariant with Da in the 
regime of small convective flow velocity. 

In Fig. 3 the maximum flame temperature (T max ) 
is plotted against the Damkohler number (Da). For 
e = 0, the T max : in gas phase is increasing with Da 


but the magnitude of temperature rise is very slight 
compared to the increment of Da. It indicates that 
the larger Da is, the more kinetical ly-strong the 
flame becomes, even though physically the fuel supply 
rate is smaller. Figure 4 is temperature distribution 
along the fuel surface which is directly related to 
the fuel evaporation rate (see Eq. (14)). As before, 
the wall temperature profiles are found to be indis- 
tinguishable from one another for a wide range of Da 
under the condition t = 0. From the above descrip- 
tion, the essential characteristics of the flame with- 
out radiant heat loss in the slow convective flow 
regime are similar to those obtained in Ref. 1. That 
is, envelope diffusion flames exist at large Da 
where the computed fuel evaporation rate, the flame 
stand-off distance and velocity profiles show certain 
similitude. 

However, the similarity relationship breaks down 
and the flame is quenched at a large Da (7.20x10®) 
when the condition e = 1, corresponding to the black- 
body surface radiation, is applied. The isotherm dis- 
tributions of the flame as a function of flow velocity 
(u„) are illustrated in Figs. 2(b) to (h). The tran- 
sition from a strong flame to extinction consists of 
a series of flame contractions. The flame fronts are 
found to maintain their positions ahead of the leading 
edge of the fuel plate and the flame quench trend by 
radiation is proceeding in the opposed flow direction 
during the transition. Comparing Fig. 2(a) (without 
surface radiation) with Figs. 2(b) and (c) we see 
that the flames shift toward the fuel plate in the 
y-direction but show little variations in x-direction. 
As u„ becomes smaller, as in the Fig. 2(d), the 
flame appears to be shortening in the x-direction. We 
start to see the complete isotherms T = 5 and T = 6 
in the gas phase and the isotherm T = 4 intends to 
form the closing end. In other words, the fire plume, 
where the fuel surface stops pyrolyzing, appears in 
the downstream gas phase. The flame size keeps 
shrinking as u„ continues to decrease. Now, we 
can see the close-ended isotherms T = 4 in Fig. 2(e) 
and T = 3 in Fig. 2(f). Finally, the isotherm T = 6 
disappears in the vicinity of leading edge zone 
(flame stabilization zone) in Fig. 2(h). And from 
Case 3s, it becomes clear that the flame contraction 
is more significant in streamwise direction than that 
in cross-stream direction. For instance, the dis- 
placement in x-direction for the flame tip of isotherm 
T = 5 is about 40 thermal lengths from Figs. 2(d) to 
2(h), but it is less than 10 thermal distances in 
y-direction. As soon as u„ reaches a critical 
value (= 1.25 cm/s), the flame no longer exists; 
radiant flame extinguishment occurs. 

Another characteristic of the transition to 
extinction can also be seen in Figs. 3 and 4. 

Figure 3 shows T max decreasing instead of increas- 
ing as u„ decreases (or Da increases) when the 
surface radiant heat loss is imposed. Near the 
extinction limit, the rate of flame temperature decay 
becomes even more dramatic. From the physical point 
of view, we have a low-power flame when the fuel 
supply rate is reduced. Figure 4 shows the tempera- 
ture profiles along the fuel surface which are coupled 
to the thermal structures in gas phase. For the con- 
venience of the study, the end of the fuel pyrolysis 
is designated as soon as Vw (blowing velocity) is 
less than 0.001 where the correspondent wall tempera- 
ture is 1.70. Since the flame front is stabilized at 
the leading edge during the transition (see Fig. 2), 
the maximum surface temperature and blowing velocity 
are expected there. In the cases of Is and 2s, 
because the flames cover the full range of the fuel 
plate x £ 50, the downstream temperature profiles 
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behind the leading edge are more or less similar to 
the ones without radiation (e = 0) but with lower 
magnitudes. From Case 3s, the fire plume starts to 
appear within the domain of interest. As u«, becomes 
smaller, the pyrolysis length becomes shorter which 
implies that the fuel supply rate is lower and the 
flame is weaker. 

Figures 2 to 4 together provide the complete 
picture for radiant flame extinguishment. The author 
mentioned earlier that the flame structure depends on 
the interaction between Da and S. From the impli- 
cation of similarity in Figs. 2(a) and 4, increasing 
Da does not affect the fuel evaporation rate as long 
as the Damkohler number is greater than 2.81x10 . 

Even the maximum flame temperature rises very little 
from Da = 2.81x10® to Da = 2.71x10® as shown in 
Fig. 3. It indicates the heat generation rate from 
combustion approaches an asymptotic value at large 
Da. As S increases with Da, the heat loss rate to 
the fuel plate becomes proportionally greater in sig- 
nificance than that of heat generation. In other 
words, at the regime of slow convective flow, the 
radiant heat loss from the fuel surface becomes so 
great that the available energy from gas phase to 
pyrolyze the fuel is greatly reduced. The fuel supply 
rate is thus reduced. As a consequence, the gas phase 
heat release rate is lowered, thereby the flame tem- 
perature is decreased, the flame size is contracted 
and the pyrolysis length is shortened. 

The large Da contributes to stabilizing the 
flame front ahead of the leading edge by providing the 
opportunity for fuel vapor to diffuse upstream into 
the flame. Since the flame front consumes a large 
portion of the limited amount of fuel vapor due to the 
shortened pyrolysis length, the extension of combus- 
tion to the downstream is delimited and the finite 
rate kinetics effect becomes crucial. Therefore the 
flame is quenched downstream and the flame contraction 
moves upstream during the transition. Finally, the 
flame disappears in the stabilization zone. 

Although the surface radiation is emphasized as 
the dominant factor for flame extinction in this 
study, its influence is coupled with the finite-rate 
reaction in the gas phase where the streamwise heat 
conduction and mass diffusion are responsible for 
flame stabilization around the leading edge and the 
downstream existence of the fire plume. Therefore 
the transition to extinction exhibits multi- 
dimensional behavior. 

Next, Case 4s is chosen to study the flame struc- 
ture within the transition. The temperature contour 
has been shown in Fig. 2(e). The flame front is sit- 
uated ahead of the leading edge of the fuel plate. 
However, the flame does not continue downstream where 
a fire plume region exists. This can be seen clearly 
in Fig. 5 which is the fuel reactivity contour distri- 
bution (-ip; see Eq. (7) for definition). The most 
reactive zone occurs in the flame front such that it 
can stabilize the flame in front of the leading edge. 
In the downstream region, the reaction cannot be 
sustained but is frozen in the gas phase due to the 
lack of fuel supply. Figure 6 gives the mass fraction 
distributions for fuel and oxidizer for this case. 
Since the most vigorous reaction zone is located at 
the flame front, a strong mass diffusion is required 
in streamwise direction to drive the fuel vapor 
upstream in order to maintain the reaction ahead of 
the leading edge. The short pyrolysis length (only 
about 10 thermal distances, see fig. 4) results in a 
reduced fuel supply which is indicated by the attach- 
ment of Yp = 0.01 fuel contour to the fuel plate. 

The oxidizer contour is shifted toward the fuel plate 
much further downstream. 


Figures 7 and 8 are relative isobar contours and 
the velocity vector field respectively. There are two 
pressure plateaus existing in the flow field. One is 
at the flame front, caused by the hot product gas 
thermal expansion, which retards/deflects the flow 
ahead of it and accelerates the flow behind it. The 
other pressure plateau is due to the existence of the 
stagnation point. With the dual effects of the 
stagnant pressure and the shear stress induced by the 
vertical blowing velocity, the flow is strongly 
deflected near the leading edge. Far from the fuel 
plate, the flow is not disturbed since the flame is 
small. In the downstream region, the reaction is 
complete and the flow is parallel to the fuel plate. 

The author is not yet aware of any experimental 
work or theoretical investigation which has the same 
model and flow configuration as the present study. 
Therefore comparisons with other works are based on 
the limited similarities which can be shown. 

The theoretical study by Sibulkin et al. (60 has 
the same solid fuel geometry as the present one but 
the flow is naturally-induced and is analyzed by the 
use of local similarity approximation. In that paper 
the authors show that surface radiation heat loss can 
result in extinction. However, the radiation is 
treated without any physical connection to the other 
factors as in the present study. Also by neglecting 
the streamwise effects due to the application of local 
similarity, the flame extinction is a local phenomenon 
which is quite different from the basic emphasis of 
this study that the flame extinguishment is a multi- 
dimensional phenomenon. 

That flame temperature is reduced because of the 
inclusion of radiation is universally recognized in 
these works by Bonne (4), T'ien (j>), Sibulkin et al . 
(6), Sohrab et al . (7_) as well as Sohrab and Williams 
[E) . The reactivity distribution in fire plume region 
resembles the profiles described in Chen and T'iens 1 
work (10) which studies the effects of the finite 
rate cFemical kinetics in the fire plume. The retar- 
dation/deflection of flow in front of the flame and 
the flow acceleration behind the flame are observed 
experimentally by Hirano and Kinoshita (11). All of 
the Navier-Stokes calculation allowing tHe thermal 
expansion, such as Lee and T'ien (_12), Mao et al. 

(13) and Ref. 1, reproduce the phenomena. Although 
tTTe flow configuration is in stagnation-point geome- 
try, T'ien (5) confirms a radiant extinguishment 
branch at small flame stretch rates where in this 
study the flame is quenched in the regime of slow 
convective flow velocity. 

CONCLUSION 

A model of a diffusion flame at the leading edge 
of a fuel plate in slow, forced convective flow under 
zero-g environment has been developed and solved 
numerically to study the flame extinguishment phe- 
nomena. The extinction mechanism in the small veloc- 
ity regime is different from that of blowoff, and 
radiative heat loss is the dominant factor. The 
mathematical system consists of two-dimensional 
Navier-Stokes momentum, energy and species equations. 

A one-step overall chemical reaction with second-order 
finite-rate Arrhenius kinetics is also included. 
Radiant heat loss on the fuel surface is applied to 
the combustion model as well as an application of an 
Arrhenius-type surface pyrolysis law. 

After the nondimensionalization procedure, two 
essential parameters are specified. These are the 
Damkohler number (Da) and the surface radiant heat 
loss parameter (S) respectively. For given fuel and 
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ambient conditions, the effect of a variable convec- 
tive flow velocity, which changes Da and S simul- 
taneously, is investigated. 

An extinction limit is identified at the small 
velocity regime. From existent envelope flame to - 
flame extinguishment, the flame front is found to 
maintain itself ahead of the leading edge of the fuel 
plate. The transition, which is a multidimensional 
phenomenon, consists of a series of flame contractions 
and pyrolysis length reductions in the opposed flow 
direction as well as a flame temperature depression. 

As the convective flow velocity is lowered to a limit- 
ing value, the radiant heat loss rate sufficiently 
outweighs the heat generation rate^and the radiant 
flame extinguishment occurs. 
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TABLE I. - NONDIMENSIONAL PARAMETERS 


Symbol 

Parameter group 

Value 

Pr 

S/S 

1 

Le 

Si D 

1 

r 

T*/T 

00 

4.33 

8 

o 

>- 


0.231 

f 


4.0 

Q 


116.8 

L 

L/CpT^ 

4.33 

E 

E/RT 

54.20 

Da 

(k*/c p uf)B 

vari able 

As 

As/p* 

2. 05x1 0 5 

Es 

Es/RT 

50.40 

S 

3 " 

eoT / p*C D U 

cd r 00 

vari able 
0, no radiation 

e 


1, blackbody 
radi ation 


TABLE II. - PARAMETRIC STUDIES 


Case 

Uoo ) 

( cm/ s ) 

Da 

S 

Is 

20 

2.81xl0 6 

0.02091 

2s 

10 

1.13xl0 7 

0.04182 

3s 

5 

4.50x10' 

0.08364 

4s 

4 

7.03xl0 7 

0.10455 

5s 

3.33 


0.12546 

6s 

2.50 


0.16728 

7s 

1.67 

4.05x10° 

0.25092 

8s 

1.25 

7.20x10® 

0.33456 (quench) 

















Figure 1.- Schematic configuration of the stabilization 

ZONE FOR AN ENVELOPE DIFFUSION FLAME (ADOPTED FROM 
REF. 1). 





cntn-c* 






Figure 2.- Concluded. 




Da 


Figure 3.- Maximum flame temperature (T max ) versus Damkohler 

NUMBER (D g ). 




Figure 4.- Temperature distr 
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